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In this paper we have compiled two new open cluster catalogues. In the first 
one, there are 119 objects with ages, distances and metallicities available, while in 
the second one, 144 objects have both absolute proper motion and radial velocity 
iy-j ! data, of which 45 clusters also with metallicity data available. 

Taking advantages of the large number of objects included in our sample, we 
CN '. present an iron radial gradient of about — 0.063±0.008 dex kpc^ 1 from the first 

sample, which is quite consistent with the most recent determination of oxygen 
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c$ ', data, and also consistent with the inside-out galactic disk formation scenarios. 

j> \ Based on the cluster sample, we also discussed the metallicity distribution, clus- 

ter kinematics and space distribution. A disk age-metallicity relation could be 
implied from those properties, although we could not give conclusive result from 
the age metallicity diagram based on the current sample. More observations are 
needed for metal poor clusters. From the second catalogue, we have calculated 
the velocity components in cylindrical coordinates with respect to the GSR for 
144 open clusters. The velocity dispersion of the older clusters are larger than 
that of young clusters, but they are all much smaller than that of the Galactic 
thick disk stars. 



gradient by nebulae and young stars, which is about —0.07 dex kpc -1 . By 
dividing clusters into age groups, we show that iron gradient was steeper in the 
past, which is consistent with the recent result from Galactic planetary nebulae 



Subject headings: open clusters: metallicity and kinematics — Galaxy: formation 
— Galaxy: evolution 



- 2 - 



1. INTRODUCTION 

Since the seminal work of Eggen, Lynden-Bell and Sandage (1962), great progress has 
been made in understanding the formation and evolution of the Milky Way galaxy. The 
progress comes, on one side, from observations concerning chemical abundances in stars 
(and clusters), gas clouds and, on the other side, from improved knowledge relevant to 
galaxy formation and evolution. 

However, some important quantities related to the chemical evolution of our Galaxy, 
such as star formation history, initial mass function, gas flow, etc., are not yet well under- 
stood. Observational data from the Milky Way disk and halo has put strong constraints on 
our understanding of those quantities. Among a variety of observables, radial abundance 
gradients along the Galactic disk is one of the most important constraints on the Galactic 
chemical evolution model. The existence of such gradients is now well established, through 
radio and optical observations of HII regions, disk stars, planetary nebulae (Henry and 
Worthey 1999; Hou et al. 2000; Chiappini et al. 2001; Maciel et al. 2002) and open clusters 
(Friel 1995, 1999). An average gradient of ~ —0.06 dex kpc" 1 is observed in the Milky Way 
disk for most of the elements, e.g. O, S, Ne, Ar and Fe. This magnitude of the observed 
gradients constrains the various parameters in the chemical evolution model, such as the 
timescales of star formation and infall (Prantzos and Aubert 1995) or any variations of the 
stellar initial mass function properties with metallicities (Chiappini et al. 2001). 

In the last decade, a number of successful models have been developed related to the 
chemical evolution of the Milky Way galaxy, but some important differences exist. One of 
them concerns the history of abundance gradient along the Galactic disk: were they steeper 
or flatter in the past? Different predictions were made by various models, although most 
of them claim that they could reproduce the majority of the observational properties both 
in the solar neighborhood and on the whole disk. Time flattening evolution is suggested by 
the models of Prantzos and Aubert (1995); Molla et al. (1997); Allen et al. (1998); Boissier 
and Prantzos (1999); Chang et al. (1999); Hou et al. (2000); Chang et al. (2002), while the 
opposite is supported by models of Tosi (1988); Samland et al. (1997) and Chiappini et al. 
(1997, 2001). 

The situation is neither settled observationally. Estimated ages of various types (PNI, 
PNII, PNIII) of planetary nebulae (PN) span a large fraction of the age of the Galaxy. 
Observations of the abundances of those objects across the Milky Way disk could, in principle, 
provide some information on the time evolution of the abundance gradient (Maciel and 
Koppen 1994; Maciel and Quireza 1999; Maciel and Costa 2002; Maciel et al. 2002). In a 
recent work, Hou et al. (2000) have made a detailed analysis for the O, Ne, S and Ar gradient 
based on the PN data of Maciel and Quireza (1999). It was shown that there is fairly good 
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agreement between model predictions and observations concerning all the properties of the 
observed abundance profiles (absolute values, gradient, scatter) for O, S, Ne and Ar. The 
model suggests that abundance gradients are steeper in the earlier epoch. However, the large 
scatter in the adopted data does not allow one to conclude on the temporal variation of the 
gradients. Nevertheless PNs suffer from large uncertainties concerning their progenitor's 
masses and lifetimes as well as their distances from Galactic center. 

On the other hand, open clusters(OCs) have long been used to trace the structure and 
evolution of the Galactic disk (Friel 1995). Since open clusters could be relatively accurately 
dated and we can see them to large distance, their [Fe/H] values serve an excellent tracer 
to the abundance gradient along the Galactic disk as well as many other important disk 
properties, such as Age-Metallicity Relation(AMR), abundance gradient evolution, disk age 
and so on (Carraro et al. 1998). 

At this point, one might ask whether the field disk populations are also able to trace the 
disk evolution. Indeed, the extensive studies by Edvardsson et al. (1993), and recently by 
Chen et al. (2000), who concentrate on disk F, G stars, show an overall radial gradient that 
is nearly independent of age. Those results are based on stars mainly restricted in the solar 
neighborhood. A more detailed analysis for the disk iron gradient was given by Cui et al. 
(2000) on the basis of 1302 field star with high resolution proper motion and parallax data 
from Hipparcos satellite. They have derived an radial iron gradient of —0.057 dex kpc -1 
within galactocentric distance from 8.5 kpc to 17 kpc. However, it is still difficult to reveal 
any pronounced gradient evolution from those results. Moreover, results from those studies 
are strongly affected by selection effects and rely on the techniques for determining individual 
stellar distances (which are heavily dependent on the adopted Galaxy potential model) that 
are much less reliable than those used to obtain cluster distances. In a recent work, Corder 
and Twarog (2001) have modelled the effects of the orbital diffusion of stars and clusters 
on the Galactic abundance gradient. The general conclusion is that the effect of diffusion 
makes a gradient shallower over time, and the cluster population offers a more viable means 
for finding detailed structure within the recent Galactic abundance gradient. 

Here, we also pointed out that our recent treatment on deriving the abundance gradient 
from open clusters in Hou et al. (2002) is in fact not proper. In that paper, we have simply 
taken four Catalogs from literatures (Carraro et al. 1998; Twarog Ashman and Anthony- 
Twarog 1997; Piatti et al. 1995; Friel 1995), and merge them just by making cross checking for 
the common clusters, without examining individually to see if there are important difference 
among clusters in the different catalogs (Twarog 2002). We refer this paper to act as a 
substitution to our previous one. 

In this paper, we compiled a set of new open cluster catalogues. The catalogue was 
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divided into two parts: CAT 1 and CAT 2. In CAT 1, we list 119 clusters with iron 
abundance, age, distance and reddening data available. This could provide statistically more 
significant information to the Galactic disk formation and evolution, such as Age-Metallicity 
Relation, abundance gradient and its time or/and spatial evolution and so on. The second 
sample consists of 144 clusters with three dimensional kinematics information available. 
From this sample, we are able to explore some statistical relations among kinematics and 
other observables. 

The paper is organized as the following: firstly in Section 2, we describe the main 
characteristic of the two samples. Then, in Section 3, we give some statistical analysis for 
the sample, mainly some metallicity and kinematics distributions. The abundance gradient 
is given in Section 4. We show that, based on our open cluster data, the abundance gradient 
of the Galactic disk was steeper in the past. In Section 5, we make some detailed discussions 
about the Age-Metallicity Relation (AMR) of the disk. Finally, a brief summary is given in 
Section 6. 



2. THE CATALOGUE 

During the past decades, a number of authors have presented their statistical studies on 
the Galactic disk based on their own open cluster catalogue. However, most of the catalogues 
suffer from either lacking of homogeneity in the cluster age and metallicity or insufficient 
three dimensional kinematic data. 

With the full release of Hipparcos Catalogue (ESA 1997) and the latest Tycho 2 cata- 
logue (H0g et al. 2000), we have seen a large growth of proper motion data for open clusters 
(e.g. Baumgardt et al. (2000); Dias et al. (2001)). A most recent compilation was given by 
Dias et al. (2002). In their catalogue, information of 1537 open clusters were present, from 
which 9% have both mean proper motion and radial velocity data simultaneously; 37% have 
distance, E(B — V) and age determinations, including 96 clusters also have iron abundances 
data available. 

We have compiled two new catalogues of the Galactic open clusters. The first one 
(hereafter CAT 1) lists 119 (including Berk29) clusters parameters for distance, age and 
metallicity. The age, distance and reddening information are all (except NGC1348, NGC2158 
and Tombaugh 2) from Dias et al. (2002), while most iron abundance data (96 clusters) were 
taken from Dias et al. (2002). The metallicities of another 23 clusters are from other 10 
literatures (Cameron 1985; Kubiak et al. 1992; Friel and Janes 1993; Friel 1995; Edvardsson 
et al. 1995; Piatti et al. 1995; Brown et al. 1996; Gratton 2000; Ann et al. 2002; Carraro et 
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al. 2002). 

Thus far CAT 1 provides a most complete open cluster sample concerning the iron 
abundance, distance and age parameters together. This sample could provide statistically 
more significant information concerning the Galactic AMR, radial iron gradient as well as 
its evolution, etc. 

In the second catalogue (hereafter CAT 2), we have listed observed kinematical data 
from literature for 144 clusters, with both radial velocity and mean proper motion available. 
The mean radial velocity data are mostly (122 of 144 objects) from a compilation in WEBDA 
database (http://obswww.unige.ch/webda/meanvr.html), primarily based on the work of 
Rastorguev et al. (1999). The absolute proper motion of 125 clusters, based on the Hipparcos 
system, are from Baumgardt et al. (2000). Mean proper motions of additional 16 cluster 
were added from compilation of Dias et al. (2001), with cluster membership probability 
derived by Tycho2 proper motions. Data of NGC2355 comes from Soubiran et al. (2000), 
and data of Coma Ber and Pleiade clusters are from Robichon et al. (1999). In fact, the above 
observed kinematic information constitutes a sub-catalogue of that of Dias et al. (2002). But, 
here in CAT 2 we have further calculated the three dimensional velocity of open clusters by 
combing with radial velocity and mean absolute proper motion data and give, for each cluster, 
the velocity components (II, 0, W) in cylindrical coordinates with respect to the Galactic 
Standard Reference (GSR). In addition, for each cluster, age and iron abundance data are 
also listed whenever available. (Notice that in the spatial velocity calculation, following 
parameters are adopted for the Sun: galactocentric distance 8.5 kpc, velocity components 
relative to LSR (10.0, 15.0, 8.0 km s _1 ) and the rotation velocity : 225.0 km s _1 . 

Our catalogue files (Table 1, Table 2) are self-explanatory. CAT 1 (Tab.l) consists of 
data for 119 (118+Berk29) open clusters. For each cluster, we list its heliocentric galactic 
coordinates in B1950.0 and the following parameters, when available: galactocentric distance; 
distance from the Sun; color excess E(B — V), age in Gyr; mean metallicity, and two reference 
codes separated by a coma, where the first number is for age, distance and E(B — V) data 
while the second number for iron abundance value. In CAT 2 (tab. 2), we provide information 
of spatial motion for 144 clusters. For each object, the following data are listed: galactic 
coordinates in B1950.0; mean radial velocity ; mean proper motions; velocity components in 
cylindrical coordinates with respect to the GSR; mean spatial velocity and error; age and 
iron abundance when available. 

In order to check if the data in our CAT 1 has any significant systematic difference 
with other published catalogues, we have made a comparison with Friel's catalogue (Friel 
1995), with 41 clusters in common. We found that the average difference in metallicity is 
less than O.lOdex, well within the typical observational uncertainty. The average difference 
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for R GC is about 0.5kpc. Note that the age indicator in Friel's work is based on the MAI 
(Morphological Age Index), which was only intended to provide a relative age ranking of 
clusters, therefore, it is not fully comparable. But there still has good overall correlation 
between Friel's catalogue and ours. 

Our following analysis will be mainly based on those two catalogues, but excluding 
cluster Berkeley 29. In CAT 1, Berkeley 29 has the galactocentric radius of 23 kpc, the 
E(B — V) of 0.15 and a metallicity of -0.18 dex, from the compilation of Dias et al. (2002). 
However, different values for these parameters of Berkeley 29 were published in the literature. 
Kaluzny (1994) gave a much smaller galactocentric distance of about 19 kpc, the reddening 
EiB — V) larger than 0.21 and based on the CMD morphology and comparison to other old 
clusters, he also deduced a [Fe/H] value of lower than -1. In the work of Noriega-Mendoza 
and Ruelas-Mayorgo (1997), who applied a new technique for simultaneous determination 
of [Fe/H] and E(B - V), a [Fe/H]=-0.30 and a E(B - V) = 0.01 were given for Berkeley 
29. As the properties of Berkeley 29 are quite uncertain we do not include this object in the 
following calculations. In general, the uncertainty for the metallicity determinations in open 
clusters is about O.ldex. 



3. STATISTICAL PROPERTIES 

3.1. The Galactic distribution of open clusters 

Using data from Dias et al. (2002) for 571 open clusters with distance and age data, we 
plotted the cluster positions on an (X,Y) coordinate system, with the zero point in X at the 
galactic center (the Sun is assumed to be at 8.5 kpc) as Figure 1 shows. Here the full line arc 
represents the solar circle about the galactic center. One will find from this figure that in the 
galactic plane, young clusters (with ages younger than that of Hyades, 0.8 Gyr, see Phelps 
et al. (1994)) distributed quite uniformly around the Sun, while roughly only 20% of the old 
clusters are inside the solar circle, most of the old ones are located further away from the 
galactic center than the Sun. This result is quite consistent with the early comprehensive 
study of Phelps et al. (1994). The deficiency of older clusters in the inner part of the disk 
has been ascribed to the preferential destruction of the clusters when they encountered with 
giant molecular clouds, which were primarily found in the inner Galaxy. 

The distribution of either old or young clusters perpendicular to the galactic plane could 
be fitted by a simple exponential law, which are plotted in Figure 2. The younger clusters are 
distributed on the galactic plane almost symmetrically about the Sun, with a perpendicular 
scale height of approximately 57 pc. In contrast, about 80% old clusters are in the outer 
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disk, outside R gc =10 kpc; this population has a scale height of about 354 pc. The derived 
scale heights are in excellent agreement with the early results of Janes et al. (1988), Janes 
and Phelps (1994). Janes and Phelps (1994) divided the open clusters into young and old 
components according to cluster's MAI, and derived their scale heights as 55pc and 375pc, 
respectively. The remarkable agreement shows that the MAI could really be a good age 
indicator for open clusters. 

The main advantage of our CAT 2 is that we have both radial velocity and mean proper 
motion available for 144 clusters. This provides a chance to probe the velocity projection on 
the Galactic plane, as shown in Figure 3. Obviously, most of the clusters are located in the 
solar neighborhood and have the velocity vectors well follow the Galactic rotational pattern. 

3.2. The metallicity distribution 

The metallicity distribution of 118 open clusters is plotted in the upper panel of Figure 
4. Here the iron abundance of about 3/4 of the OC sample (with [Fe/H] > -0.2) has roughly 
a Gaussian distribution which peaks at the solar value. Meanwhile a metal poor metallicity 
tail is also clearly seen. Here we divide our OC sample by two groups, that is, Metal-Poor 
(hereafter MP) component and Metal-Rich (hereafter MR) component, with the dividing 
line between them, somewhat arbitrarily, at [Fe/H]= —0.2. In the lower panel of Figure 4, 
we show two histograms for the open clusters of the above two groups. By assuming an 
exponential law, we can derive their scale-heights to be 535 pc and 106 pc, for the MP and 
MR components, respectively. Taking the scale heights of the Galactic thick and thin disk 
as 760 pc and 260 pc (Ojha et al. 1996), we could see that spatially, the MP group might 
be either within the tail of the thick disk or in the outskirts of the thin disk while the MR 
group is just a thin disk component. However, our OC metallicity sample is surely not a 
complete one, and is subject to a variety of observation effects. For example, the outer disk 
clusters are subject to significant selection effects - they can be seen more readily if they are 
at larger distances from the plane, and they are likely to be older, since younger clusters 
which live closer to the plane will not be as visible, and more difficult to observe. 

This can also be seen clearly in Fig. 5, where we have plotted the dependence of clusters 
vertical height on the galactocentric distance for the whole sample with distances data avail- 
able (Fig. 5a) and for those also have metallicity data (Fig.5b). In Fig. 5b, we see that most 
of the MP OCs are outer disk objects (with a median R gc ~ 11.3 kpc and the majority(65%) 
are older than 0.8 Gyr, see also Fig. 8,9) and relatively far away from the Galactic plane, 
with a median distance about 326 pc. For MR clusters , most of them are in the inner disk 
(Rgc <10 kpc), distributing in the immediate solar neighborhood (with a median R gc ~ 8.7 
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kpc), and about 78% are young (with ages < 0.8 Gyr) objects (also see Fig. 8,9). This 
phenomena is very likely to imply the possible existence of the age metallicity relation in 
open clusters. 

Meanwhile, when R gc <10 kpc, most OCs, no matter MP or MR, are much closer to the 
Galactic plane, with a median height z ~ 84 pc. Especially, there is few clusters observed 
in the region of Rgc <6 kpc. This has been attributed to the destructive power of the large 
numbers of giant molecular clouds in the inner regions of the Galaxy (van der Bergh and 
MMcClure 1980). There is much evidence that leads us to believe that open clusters have 
been selectively destroyed near the plane of the disk and only those clusters whose orbits 
keep them away from the Galactic plane can survive long enough to appear as outer disk, or 
metal-poor clusters. In the meanwhile, this also imply that part of the outer clusters (they 
are metal-poor and with high-z) might be formed as a result of disturbances to the Galactic 
disk, possibly caused by tidal interactions with other galaxies or infalling gas, as suggested 
by Janes and Phelps (1994). However, it is still hard to understand why there are almost no 
high z clusters in the region of Rgc ~ 6.5-8 kpc, compared with the outer disk results. 

As a comparison, both globular and open cluster metallicity distribution are plotted 
in the upper and lower panels of Figure 6. We can see a clear overlap between metal rich 
globular cluster and metal poor open clusters around [Fe/H]~ —0.4 dex. If the age metallicity 
relation do exist, then this could be another evidence which support the idea that a possible 
connection between the halo and disk population exist, both in their chemical and dynamical 
history. 

Based on a survey of proper motion stars, Carney et al. (1990) pointed out that the 
Galactic halo population had the chemical and dynamical history almost independent of the 
disk. However, from a study of the oldest open clusters, Phelps et al. (1994) found that the 
oldest open clusters (Bel7 with an age about 12.5Gyr) have the ages compatible with that 
of the youngest globular clusters, suggesting that there may have been little or no delay in 
time between formation of the halo and the onset of the development of the disk. Phelps's 
argument rests largely on the age of Bel7 clusters, however, recent works have given an age 
of about 9 Gyr to this clusters (Carraro et al. 1999), and so there does still appear to be a 
gap between the formation of the halo and the thick disk. On the other hand, if the cluster 
metallicity is related to their age, the metallicity overlapping of the clusters, as we have 
presented, may be another indication of the connection proposed by Phelps et al. (1994). In 
fact, the distinction between "open" and "globular" (so-called "super") clusters may turn 
out to be largely an artificial one (Larsen 2002). They could be both originated from the 
Super Star Clusters (SSCs) which are observed in large numbers in interacting galaxies and 
merger remnants. Our Milky Way disk is very likely undergone a process of minor mergers 
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in the early epoch. The thick disk is plausibly the result of heating of the thin disk through 
such events (Wyse 2001). 

3.3. Metallicity vs. kinematic 

In Figure 7, we plot the relation between cluster rotational velocities around the Galac- 
tic center, 6, and the galactocentric distance. There exhibits,from least-square fitting, an 
insignificant slope of about -2.5 km s -1 kpc -1 , with quite large scatter. In the right panel, 
we present the dependence of velocity dispersion on the cluster age. The young clusters have 
a smaller velocity dispersion, which is expected from their small z scale height. Although 
the velocity dispersion for the older clusters is about 20 km s -1 , it is still much smaller 
than that of the Galactic thick disk stars, which is about 50 km s -1 . Anyway, for most of 
the CAT. 2 clusters ( 90%), their heights from the galactic plane are well within 200pc, they 
are just thin disk objects. 

In our kinematic sample, there are only 2 objects in the outer disk (Berkeley 31 and 
Dolidze 25, both with R gc >14 kpc). They have unreliable proper motion and radial velocity 
results (with relative errors up to about 50%) and thus were not included in the above radial 
gradient fitting. 

4. DISK METALLICITY GRADIENT FROM OPEN CLUSTERS 

4.1. The abundance gradients 

The first radial metallicity gradient using open clusters was given by Janes (1979) based 
on DDO and UBV photometric data of 41 disk objects (part of them are field stars). The 
derived gradient is —0.05 dex kpc -1 . Panagia and Tosi (1981), by matching theoretical 
isochrones to HR diagrams of 20 clusters with age less then 1 Gyr, derived an iron abundance 
gradient of —0.095 dex kpc -1 . A similar result was also obtained by Cameron (1985) based 
on 37 clusters with mixed ages. By introducing a weighting system in order to evaluate and 
compare the published parameters in Lynga (1987) Catalog of Open Clusters data, Janes et 
al. (1988) determined some basic parameters of 413 open clusters, such as ages, distances, 
linear diameters and so on. Among them, 87 open clusters have the metallicity data. They 
have derived a gradient about —0.133 dex kpc -1 . By separating the clusters into age groups, 
they found that young clusters has much smaller gradient than that of older clusters. Besides, 
all those authors found some indications that the gradient became shallower in the direction 
of the Galactic center and steepen in the outer parts of the Galaxy. 
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Friel and Janes (1993), (hereafter FJ93) presented their results from a spectroscopic 
study of a sample of giant stars in 24 open clusters. They derived a galactocentric radial 
abundance gradient of [Fe/H] about —0.088 dex kpc -1 . A subsequent revision of the FJ93 
result was presented by Friel (1995), using additional spectroscopic results and a more uni- 
form set of cluster properties. From a sample of 44 clusters, Friel (1995) derived an iron 
gradient of —0.091 dex kpc -1 . At the same time, Piatti et al. (1995) derived a much smaller 
gradients, —0.07 dex kpc -1 , from a sample of 63 open clusters with a wide range of ages. 
These results are quite consistent with the recent result of Friel (1999) who obtained a gradi- 
ent about —0.06 dex kpc -1 . Another gradient result was presented by Carraro et al. (1998) 
recently. The metallicities of all selected 37 clusters were obtained spectroscopically. The 
final gradient was about —0.085 dex kpc -1 , agreed with earlier result of FJ93. By dividing 
the sample into age bins, it was found that the present-day gradient is a little shallower than 
the past one, while the middle epoch seems to display a steepening of the gradient. 

The presence of a linear gradient for open clusters has been questioned by Twarog 
Ashman and Anthony- Twarog (1997) (hereafter TAA97). TAA97 put forth an alternative 
description, namely, step function, about the radial abundances distribution of the open 
clusters. Within this work, a set of 76 clusters with abundances based upon DDO and/or 
moderate dispersion spectroscopy has been transformed to a common metallicity scale and 
used to study the local structure and evolution of the Galactic disk. They found that the 
metallicity distribution of clusters with galactocentric distance is best described by two 
distinct zones, with a sharp discontinuity at Rqc = 10 kpc. Between Rqc = 6.5 kpc and 
10 kpc, the clusters have a mean metallicity of 0.0 dex with, at best, weak evidence for a 
shallow gradient over this range, while those beyond 10 kpc have a mean value about —0.30 
dex. This two-step distribution seems quite similar with the nebula results of Simpson et 
al. (1995). Neglecting this two-step phenomena, a least square fitting results in a gradient 
about —0.067 dex kpc -1 between 6 and 15 kpc if cluster BE21 was excluded because both 
metallicity and distance of this object are quite uncertain. 

The existence of radial iron abundance gradients is also confirmed by our new up-to-date 
sample. The result is shown in the upper panel of Figure 8. By equal-weighted least-square 
fitting, we derived a radial abundance gradient of —0.063 ± 0.008 dex kpc -1 , which agrees 
well with most of the previous open cluster results. And it is also similar to the gradients 
obtained from other tracers, such as disk HII regions and planetary nebulae (see a summary 
in Hou et al. (2000) ). The existence of gradient along the galactic disk provides good 
opportunity to test theories of disk evolution and stellar nucleosynthesis. It suggests that 
the role of the Galactic bar in inducing large scale radial mixing and therefore flattening the 
gradient has been rather limited; alternatively, the bar could be too young(<lGyr) to have 
brought any important modifications to the gaseous and abundance profile. However, we 
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must notice that our current knowledge on the iron gradient as derived from open clusters 
is far from being clear. Open clusters span a wide range of age, from several millions years 
to several Gyr, therefore they do not trace the young component of the galactic disk. The 
result we obtained is somewhat an averaged one (over age). The obtained similarity of 
gradient between iron and other elements, such as oxygen, is quite surprising since the sites 
of nucleosynthesis for iron and oxygen are quite different. It is well know that iron is mainly 
produced in type SNIa, while oxygen is largely a product of SNII, that is from massive stars. 
So the abundance history is very different for those two types of elements. The gradient 
similarity might be simply coincidental, or further investigations should be given in the 
production nature of those elements. 

We also derived a vertical abundance gradient of —0.295 ± 0.050 dex kpc -1 (lower panel 
of Figure 8). This is consistent with the result of Carraro et al. (1998). 

4.2. Gradient evolution in the Galactic disk 

As we have pointed out in Section 1, the behavior of gradient evolution along the Galactic 
disk is a major problem for different chemical evolution models. Open cluster system is an 
ideal template for this analysis because OCs have relatively well determined ages, distances 
and metallicities. 

In the upper panel of Figure 9, we show gradients for two sub-samples with cluster age 
< 0.8 Gyr (80 clusters) and > 0.8 Gyr (38 clusters), respectively. The fitting results are 
—0.024 ± 0.012 dex kpc -1 for younger clusters, —0.075 ± 0.013 dex kpc -1 for older ones. If 
we take the mean age for the youngest and oldest clusters as 0.00 Gyr and 6.00 Gyr (this is 
somewhat arbitrary, just for illustration purpose) in our sample, we can estimate an average 
flattening rate of 0.008 dex kpc~ 1 Gyr~ 1 during the past 6 Gyr. Similar value is obtained 
by Maciel et al. (2002) from PN data for [O/H]. 

As we have indicated in the Introduction, that the time evolution of the abundance 
gradient along the Galactic disk is crucial in discriminating different theoretical models that 
adopt various prescriptions used for the time dependence of the SFR and the infall. Our 
current open clusters sample could surely provide some insights on this subject. The time 
flattening tendency we obtained supports the 'inside-out' disk formation scenarios with infall 
time scale dependent on radius from the disk center (Boissier and Prantzos 1999; Hou et al. 
2000; Chang et al. 2002). 

In the lower panel of Figure 9, we divided clusters into inner ( < 10 kpc) and outer 
groups. The corresponding gradients are —0.040 ± 0.022 dex kpc -1 and —0.047 ± 0.023 
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dex kpc -1 , respectively. We can see that the inner disk exhibits roughly the same (or a 
bit smaller) gradient as the outer part. This result is also consistent with the abundance 
gradient determined by using Cepheid in the solar neighborhood (Andrievsky et al. 2002). 
However, in our CAT 1, the inner most cluster is located at a galactocentric distance about 
6.8 kpc, it is necessary to have more inner clusters data (between 3 kpc and 7 kpc ) in order 
to further check the gradient behavior for the inner disk. If the Galactic bar does play the 
role, then the inner gradient could be more flat compared with outer part. 

Our cluster sample is nearly 50 % more than that of TAA97, and we did not find evidence 
of any abrupt discontinuity. A similar conclusion was reached by Friel (1999), using high 
resolution abundance determinations for metallicity calibration. 

5. DISK AGE-METALLICITY RELATION 

The age- metallicity relation(AMR) for the Galactic disk provides useful clues about the 
chemical evolution history of the Milky Way, and also put an important constraint on the 
theoretical models of the disk. The observed abundance data generally show a decrease of the 
stellar metallicity with increasing stellar age, indicating a continuous growth of the metals in 
the ISM during the life of the Galaxy. The early study on AMR for nearby stars by Twarog 
(1980) found that the mean metallicity of the disk increased by a factor of five between 12 
and 5 billion year ago and has increased only slightly since then. This was also confirmed by 
latter photometric survey of Meusinger et al. (1991). With the high resolution spectroscopic 
data, Edvardsson et al. (1993) showed a plot of iron abundance versus relative ages for the 
189 stars in the solar neighborhood. The overall trend of a slowly increasing abundance with 
decreasing age was consistent with the previous photometric results. However, the most 
striking feature of their result is the large scatter around the average trend, which marks a 
weak correlation between age and metallicity. This spread was, as they pointed out, in part 
due to selection bias for the programme stars, and at least partly intrinsic, since the mean 
errors in [Fe/H] measurement and logarithmic age derivation are much less than the scatter. 

In a recent paper, Feltzing et al. (2001) have re-examined the Galactic AMR in the solar 
neighborhood based on a sample of 5828 dwarfs and sub-dwarfs from Hipparcos Catalogue. 
They found that the solar neighborhood age-metallicity diagram is well populated at all 
ages and especially that old, metal-rich stars do exist, which have been omitted in previous 
samples. This indicates a complete lack of enrichment over the age of Galactic disk among 
the fields stars in the solar neighborhood. 

Using open clusters to explore the AMR has the main advantage both in abundance and 
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age determinations since one is dealing with a group of stars and the result is less susceptible 
to individual errors (Carraro et al. 1998). Cameron (1985) was the first to probe the AMR 
from open cluster data, and found no age-metallicity relation based on his cluster sample. 
This is not surprising since the metallicity of the Galactic disk increased only slightly during 
the past 5 Gyr, while his sample of 38 clusters contained no objects older than 5.1 Gyr. 
More recently, Carraro et al. (1998) compiled a relatively homogenous sample of 37 open 
clusters. The data have more expanded cluster ages up to 9 Gyr. After correcting for the 
radial abundance gradient, the derived AMR showed similar trend to that of nearby stars. 

In this paper, we have present a new open cluster catalogue with much more objects. The 
results, based on this larger sample, would be statistically more reliable. As we have shown in 
the Sect 3.2, statistically, the space distributions (scale heights for metal poor and metal rich 
groups, for young and old clusters) of open clusters are very likely imply the existence of age- 
metallicity relation in the Galactic disk. In Fig. 10, we plot the dependence of metallicity on 
the cluster age, after correcting the radial metallicity gradient. Unfortunately, it is difficult 
to draw any conclusive indication for AMR based on this plot due to the deficiency of very 
older clusters. More observational efforts should be added in finding more older clusters. 

The significant spread of the AMR seems real, but its origin is not yet clear. For 
the scatter in the AMR of nearby stars, many possible causes have been suggested, such 
as orbital diffusion of stars, inhomogeneous chemical enrichment in the Galaxy evolution, 
overlapping of different galactic substructures and so on. All the above mentioned effects 
may contribute the observed scatter, while for open clusters, the result should not be very 
sensitive to orbital diffusion effects (Corder and Twarog 2001). Therefore, the scatter of 
AMR along the Galactic disk from both clusters and field disk stars is an essential feature 
in the formation and evolution of the Milky Way. 

6. Summary 

The main work of this paper is to compile a most complete open clusters sample with 
metallicity, age, distance data as well as kinematic information available. And upon this 
sample, some statistical analysis on spatial and metallicity distributions have been made. 

We derived an iron radial gradient about — 0.063±0.008 dex kpc -1 from the CAT 1, 
which is quite consistent with the most recent determination of oxygen gradient in nebulae 
and young stars. By dividing clusters into age groups, we show that iron gradient was steeper 
in the past, which is consistent with the recent result from Galactic planetary nebulae data. 
Our result supports the inside-out Galactic disk formation mechanism that invoking SFR 
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and infall time scale various with radius. 

The spatial variation of gradient was also explored. When the clusters were divided into 
inner and outer groups, we found that the abundance gradient was a bit shallower inside 
10 kpc. But the inner most cluster in our sample is in Rqc = 6.8kpc, we need more data 
of clusters in the inner region. This could be helpful to judge the radial flow effect on the 
current galactic chemical evolution model. 

From scale heights of metal poor and metal rich clusters, we noticed that the metal- 
licity could be related to the age. However, by plotting directly the dependence of cluster 
abundancess on their ages, no strike slope in AMR was found. However, the paucity of 
metallicity of very old open clusters made it impossible to give a definite conclusion based 
on the current sample. 
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Fig. 1. — Spacial distribution of the open clusters on the galactic plane. The open circles 
are old clusters with age great than that of Hyades (0.8 Gyr), and the dots are for younger 
ones. The Sun is at X = 8.5,Y = kpc. The galactic center is at (0,0). The circle has a 
radius of 8.5 kpc, centered on the galactic center. 
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Fig. 2. — Number distribution from the galactic plane, z, for old and young sub-groups of 
open clusters. The fitted scale heights for the two groups are 354pc and 57pc, respectively. 
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Fig. 3. — Velocity projection on the Galactic plane for 144 clusters, which have both radial 
velocity and mean proper motion available. The Sun is located at X = 8.5 kpc, Y = pc. 
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Fig. 4. — Upper panel: metallicity distribution of the Galactic open clusters. Lower panel: 
number distribution from the galactic plane, Z, for metal poor and metal rich sub-groups of 
open clusters. We can see that the two groups are different in their scale properties, with 
scale height about 535 pc and 106 pc,respectively. 
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Fig. 5. — Upper panel: height from the galactic plane vs. galactocentric distance, for all 571 
clusters with data available. Lower panel: same as the upper panel, but only for 118 OCs 
with also iron abundance data. It can be seen that the most metal poor clusters are located 
in the outer part of the Galactic disk. 
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Fig. 7. — Left: correlation between clusters rotational velocity 6 and galactocentric distance. 
Right: velocity dispersion vs. cluster age. 
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Fig. 8. — Radial (upper panel) and vertical (lower panel) abundance gradient for 118 open 
clusters. The least-square fitting results in a gradient of —0.063 ± 0.008 dex kpc -1 and 
—0.295 ± 0.050 dex kpc -1 , respectively. The typical error bar for [Fe/H] is about O.ldex, 
as showed in lower left corner of the figures. When deriving the vertical gradient, the radial 
gradient has been corrected. 
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Fig. 9. — Upper panel: Time evolution of the iron gradient. Triangles show clusters with 
age less than 0.8 Gyr, stars show clusters with age greater than 0.8 Gyr. The gradients are 
—0.024 dex kpc -1 and —0.075 dex kpc -1 , respectively. Lower pane: gradients for inner disk 
(within 10 kpc) and out disk clusters. The corresponding gradients are —0.040 dex kpc -1 
and —0.047 dex kpc -1 , respectively. 
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Fig. 10. — Age-Metallicity Relation (AMR) for the 118 open clusters after correcting for the 
radial gradient. The solid line is a least-square fitting for the open cluster data. 
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Table I. Age, Distance and Metallicity of Open Clusters 
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Table 1 — Continued 
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We list 119 clusters, while in the paper, Berkley 29 was not included in our calculations. 
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Table 2. Kinematic Parameters of Open Clusters 
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